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1.  INTRODUCTION 

There  is  a  naturally  conjectured  inequality  for  sums  of  random 
variables  that  seems  hard  to  prove.  Consider  Y^,  1  <_  i  <_  n,  indepen¬ 
dent  and  identically  distributed  positive  random  variables.  Let 


W 


6lV92Y2+- 


•+e  y 

n  n 


with  0^  positive  constants.  The  conjecture  is  this:  among  all  ways 

of  varying  9^  that  preserve  the  mean  (so  0j+»»»+0  stays  fixed), 

the  tails  of  W  are  smallest  when  all  0^  are  equal. 

It  is  not  hard  to  show  that  the  variance  of  W  (or  the  expected 

value  of  any  convex  function)  is  smallest  when  all  0^  are  equal.  It 

often  seems  as  if  much  more  is  true. 

For  symmetric  distributions  a  fairly  general  result  is  known. 

Let  f  be  the  common  density  of  Y. , . . . ,Y  .  Define  0  =  (0,  ,6.,  . . .  ,9 

in  '12 

Proschan  (1965)  showed  that  if  f  is  symmetric  (about  zero)  and  log  f 
is  concave,  then  P[W  t]  is  Schur-convex  in  0  for  t  >  0. 

For  positive  random  variables  (having  asymmetric  distributions) 
very  little  is  known.  We  have  chosen  to  work  mainly  with  the  gamma 
distribution.  Our  results  are  summarized  in  the  following. 

Let  Y^  have  a  Gamma  (a,  0)  density  f(y)  =  e  9^ya  ^Ba/r(a)  on 

(0,®).  For  n=2  we  show  - - 

j  Accc  T  -  t, 

r  - 

!  :r ' 

P[W  <_  t]  is  Schur-convex  in  6  if  |  ; 


t  <  a  (9. +9.)/6,  and 

—  12  i 


E.  •* 


P[w  >  t]  is  Schur-convex  in  0  if 

t  >  (a  +-^)  (G1+02>/6. 

Specializing  to  a  standard  exponential  with  a=8=l,  both  tails  are 
smallest  when  0^=r„.  Here  the  left  tail  starts  any  place  left  of 
the  mean  (=  anc*  r^SHt  tail  starts  any  place  to  the  right 

of  the  mean  times  3/2.  No  simple  convexity  holds  between  these  two 
values . 

For  general  n  >  2,  our  results  go  in  the  same  direction  although 
they  are  cruder.  For  a  sum  of  n  Gamma(u,6)  variables 

P[W  _  t ]  is  Schur-convex  for  6  in  the 

region  {0:  min  0  _>  tB/(na+l)},  and 

lli<n  1 

P[W  2.  t]  is  Schur-convex  in  ~  for 

t  >  (na+1) (r.+0o+* • *+0  )/B  . 

—  1  z  n 

Our  results  cover  positive  linear  combinations  of  chi-square  variables 

These  arise  in  the  asymptotic  distributions  of  nonparametric  goodness  of 

fit  tests.  See  Chernoff  and  Lehmann  (1954)  and  Moore  (1978)  and  Alvo, 

Cabilio  and  Fiegen  (1982)  (for  the  average  Kendall  tau  statistic).  The 

results  of  this  paper  apply  because  a  linear  combination  of  the  form 
2  2 

c^xn  +  c2^m  be  written  as 


where  0^  =  2c^  for  1=1,..., n  and  0^  =  2c2  for  1 = n+1, . . . ,n+m, 
and  the  independent  gamma  variables  each  have  shape  parameter 

a  *  and  8  =  1. 

Positive  linear  combinations  of  chi-square  variables  also  arise  when 
positive  linear  combinations  of  sample  variances  from  normal  populations 
are  formed. 

Linear  combinations  of  chi-square  variables  can  also  occur  as 
limiting  distributions  of  U-statistics .  For  example,  let  h(x,y)  be 
a  symmetric  function  with  associated  U-statistic 


U„ST  I  h(X.,X.) 


(£>  i<J 


where  the  X^  are  independent  identically  distributed  random  variables 

with  distribution  function  F.  Assume  E(h  (X^,X2>)  is  positive  and 

finite.  Define  h^(x)  to  be  E(h(x,X^))  and  to  be  the  variance 

of  h1 (X0) .  Let  0  =  Eh(X,,X„).  If  5,  is  zero,  then  n(U  -9) 
x.  t  i  l  l  n 

converges  in  distritubion  to  a  constant  plus  a  linear  combination  of 
independent  chi-square  random  variables  with  one  degree  of  freedom. 
Details  and  examples  can  be  found  in  Gregory  (1977)  or  Serfling  (1980) 
Section  5.5.2. 

Our  results  may  be  applied  to  study  linear  combinations  of  exponen¬ 
tial  variables.  These  arise,  for  example,  as  waiting  time  distributions 
in  pure  birth  processes  (see  Feller  (1971),  page  41). 

We  also  give  results  for  the  Weibull  distribution  when  n=2:  Let 
the  Y^’s  have  a  Weibull  (S,y)  density 


We  show  that  P[W  <_  tg]  is  Schur-concave  in  8  if 

1 

(a)  tQ  _>  (91+62)((1  +  YjjjO/Y)8  when  S  _>  1 

or 

i  ?  — e- 

0>)  tQ  i  (e1+^2)  ( (i  +  ^)/y)  2  s  when  0  <  e  <  i  • 

Our  results  show  that  the  tail  of  the  linear  combination  W  is 
minimized  when  all  the  coefficients  8  are  equal.  This  is  analogous 
to  results  for  binomial  random  variables  in  the  master's  thesis  of 
Chebychev  and  in  the  work  of  Hoeffding  (1956)  and  Gleser  (1975)  which 

are  described  in  Section  K  of  Chapter  12  of  Marshall  and  Olkin  (1979). 

(Our  results  resemble  Gleser's  but  the  inequalities  go  in  the  opposite 
direction.)  We  believe  similar  results  hold  much  more  generally,  e.g., 
for  all  naturally  occurring  exponential  families.  We  have  been  unable 
to  provide  a  general  theory. 

It  is  possible  to  consider  the  question  of  Schur-convexitv  or 
concavity  in  other  metrics  for  the  0^'s  and  it  should  be  noted  that 
Tong  (1980)  has  shown  for  gamma  random  variables  that  P(W  _<  t)  is 
strictly  Schur-concave  in  (9^, . . .  ,8  *")  for  all  positive  t  and 
n  >  2. 

In  Section  2  of  this  paper  we  treat  the  case  of  a  positive  linear 
combination  of  two  independent  identically  distributed  gamma  random 
variables.  In  Section  3  we  look  at  this  case  for  a  Weibull  distribution. 


In  Section  4  we  examine  the  distribution  function  of  a  positive  linear 
combination  of  n  independent  identically  distributed  gamma  random 
variables. 

2.  POSITIVE  LINEAR  COMBINATIONS  OF  TWO  GAMMA  VARIABLES 

In  this  section  we  set  n=2.  Without  loss  of  generality  we  assume 
=  1  and  assume  the  scale  factor  B  is  one  in  the  Gamma  (a, 6) 
distribution.  Then  with  p  =  6^  we  have 

W  =  pY:+  <l-p)Y2 

where  0  <  p  <  1  and  the  Y^  are  independent  random  variables  each  with 
density 

f  (x)  “X01  1  e  X/7  (a)  , 

for  x  >  0,  where  a  is  positive. 

In  the  case  n=2  the  Schur-convexity  or  Schur-concavi tv  in 
of  the  distribution  function  of  W  is  determined  by  examining  the 
sign  of  the  following  derivative 

^  P(pV  <1-P>Y2  i  t)  , 

where  j  <  p  <  1. 

We  show  for  all  p  there  exists  a  value  tg(p)  such  that  the 
derivative  is  positive  for  t  less  than  tg(p)  and  negative  for  t 


greater  than  tQ(p).  Also  tQ(p)  lies  between  a  and  a  +  y. 

Thus  P  (W  <_  t)  is  Schur-concave  in  6  for  t  <_  a  +  j  and  Schur- 
convex  for  t  £  a.  Furthermore  P(W  <_  t)  is  neither  Schur-concave 
nor  Schur-convex  in  9  for  a  <  t<a+  y.  This  result  was  observed 
empirically  by  Solomon  (1961)  in  distribution  function  tables  for  the 
case  of  a  =  y  (when  the  gamma  variable  is  a  multiple  of  a  chi-square 
variable  with  one  degree  of  freedom).  The  result  was  shown  analytically 
by  Diaconis  in  unpublished  results  for  a  =  1  (when  the  gamma  variable 
is  an  exponential  random  variable).  See  Marshall  and  Olkin  (1979), 
Chapter  12,  Theorem  K.3. 

Notation:  For  the  rest  of  this  section  Y  ,Y2  are  independent 

a-1  -x 

gamma  random  variables  with  density  x  e  /T(a)  and  p  will  always 
satisfy  <  p  <  1. 

Furthermore,  we  will  often  write 


F  (t)  =  P[pY.  +  (l-p)Y,  <  t]  . 


Lemma  1 :  The  sign  of 


jjPIpYj+U-pttjit) 


is  the  same  as  the  sign  of  S  (t)  defined  by 

P 


S  (t)  =  [  (p-x)x31-1  (1-x)  1-1e  txc 

p  Jo 


v'a'-.- vV \\ 


where  X 


In  fact 


=  ±XZ±. 

p(l-p)  * 


P[pY1+  (l-p)Y2<  t]  =  KSp<t) 


where  K  = 


.2  /+1  -t/(l-p) 


Ip  Ci-p) 


Proof : 


implies 


(tr(a))‘ 


Fp(t)  . 


p  a-1  -y 
r  y_  e 


(a) 


(1-p)  a-1  -u 

u  e 


T(a) 


tFp(° 


P  <t-y)_  yU  Vy  (t-py°  1  _e 

o  a-p)2  r(a)  r(a 


Making  the  change  of  variable  x  =  ^-v  yields 


4-  p  (t>  =  £ 

dp  p  p 


fl 

-0 


(t-tx/p) 

(1-p)2 


a-1 


( 


t-tx 

1-p 


) 


a-1 


f  t2  f*1  e~  a-P}  f1 
tp (1-p) j  (tr(a))4  j0 


(p-x) [x(l-x) ] 


where 


Remark:  The  next  lemma  shows  that  for 


p  <  1,  the  function 


e  PCSp(t)  is  strictly  decreasing  in  t  and  thus  can  have  at  most  one 

change  of  sign.  From  the  definition  of  S  (t)  it  is  clear  that 

P 

S  (0)  0  and  that  S  (t)  <  0  for  all  sufficientlv  large  t.  Thus 

P  P 

there  exists  Cq(pj  0  such  that 


r-  F  (t)  >  0  for  t  <  t,.(p)  and 

•p  p  U 


f  (t)  <  0  for  t  >  t-.(p)  . 
?  0 


1  —  \  p  £ 

Lemma  2:  For  —  <  p  <  1,  the  function  e  S^(t)  i s  strictly 


decreasing  in  t. 


Proof : 


..-'P's  (t)  -  f1  (p-,).*-1!!-.]"1,  l(x-p)dx 


so  that 


T7  ie  ptS  (t)] 


which  is  negative.  QFD 


,  /  \2  '1-1  , ,  S  1-1  '  t(x-p) 

-  X  i  (p-x)  x  (1-x)  e  dx 


The  next  lemmas  3,  A,  and  5  are  used  in  the  proof  of  Theorem  1 


which  states  that  tn(p)  is  less  than  (a  +  -5-) 


Lemma  3 : 


3p 


Vc> 


has  Che  same  sign  as 


1  -t/h(s) 

e _ 

0  (h(s))2a+1 


[  (l-2s)sa_1(l-s)Cl  1  ] ds 


where  h(s)  =  ps +  (1-p) (1-s) • 

Proof : 


a-1  -x  ct-1  -v 
x  e  y  e  -  ,  , 

~n IT  rlaT  dxdy  ’ 

where  integration  is  over  the  values  such  that  x  _>  0,  y  >_  0 
px  +  (l-p)y  <_  t . 


Fp ( t )  - 


ff 

jj 


Substituting  r  =  x+y  and  s  =  x/(x+y)  this  becomes 


F  (t)  =  r(a)~2  I  sa  1(l-s)a“1ds 

P  Jo 


t/h(s) 


2a-l  -r 
r  e  dr 


so  that 


f2j-  r1  „-t/h(s)  ,  , 

~  F  (t)  =  -£— x  — - j-t  [  (l-2s)sJ-1(l-s)^ 

3p  P  f (a) 2  -  0  (h(s))2a+1 


Lemma  A:  Define 


f (x,y) 


log(^) 


log  (^) 


X  V 


and 


]  ds  . 


for  x  ^  v  and  f(x,x)  =x.  Then  f(x,v)  is  strictly  Schur-concave 


Proot ;  Note  that 


T7 )f(x,y)  *  (r-^-)  -  [x--+y  ,!  log (y /x) 
y-x  Wv-x'r 


So  tor  x  v  we  must  show  that 


x+v 

v-x 


v(y-x)2j 


log(v/x) 


or  equivalently 


2_  2 

‘  2  %  <  log(y/x) 
x  +y 


or  equivalently  (substituting  u  =  y/x) 


—j —  <  log  u  for  u  >  1  . 
u+1 


Both  sides  are  equal  when  u=l  so  it  suffices  to  differentiate 
both  sides  and  show  that 


4u 


(u‘+l) 


2  <  Z  ' 


tor 


u  >  1 


This  is  immediate.  OED 

Lemma  5 :  For  x  >  0  define  g(x)  =  "  +  log  x. 

(a)  If  0  <■  v  <  z  and  v+z  <_  2,  then  g(y)  >  g(z). 

(b)  If  g(v-t-z)  _>  g(y-z),  then  g(v+u)  >  g(y-u)  for  0  <  u  < 


r»i«] 


(a)  For  y  <  z,  g(v)  >  g(z)  is  equivalent  to  f(y,z)  <  1  where 
f  is  the  function  from  the  previous  lemma.  By  Schur-concavity 

f(y,z)  <  f(^,^0  =  (y+z)/ 2 

so  that  g(y)  >  g(z)  when  (y+z)/2  1. 

(b)  For  r  <  s,  g(r)  <  g(s)  if  and  only  if  f(r,s)  >  1. 
Similarly  g(r)  <_  g(s)  if  and  only  if  f(r,s)  >_  1.  Thus 

g (y+z )  >_  g(y-z)  implies  f(y-z,  y+z)  _>  1.  Now  Schur-concavity  yields 
f(y-u,y+u)  >  f(v-z,y+z)  _>  1  for  0  <  u  <  z.  Thus  g(y+u)  >  g(y-u) 
as  desired.  QED 

Theorem  1:  (See  Remark  preceding  lemma  2). 

r-  1 

For  p  > 


tQ(p)  <  a  +  |  ; 

3  1 

that  is,  —  F  (t)  <  0  for  t  >  a  +  —  . 

dp  P  —  2 

Proof :  We  use  the  result  of  lemma  3.  The  function 

Ct“l  ct  — 1  1. 

(l-2s)s  (1-s)  is  odd  about  the  value  s  =  ^*  Thus,  to  show  that 

£ 

t  F  (t)  <  0,  it  suffices  to  show  that 
dp  P 

-t/h(s)  -t/h(l-s) 

— - - -  <  — - 

(h(s))2a+1  (h(l-s))2a+1 


for  0  <  s  <  2’  By  taking  logs  and  doing  some  algebra,  this  is  seen 
to  be  equivalent  to 


c 

h(s) 


+  log 


h(s)] 

.  c  J 


> 


c 

h(l-s) 


+  log 


where  c  =  • 

=  ps  +  ( 1 — n )  (1-s ) 
will  hold  if 


Now  use  lemma  5(a).  Since  the  function  h(s) 
is  increasing  for  p  >  y,  the  previous  inequality 


h(s)  +  h(l-s)  <  2 
c  c  — 


But  h(s)  +  h(l-s)  =  1  for  all  s,  so  this  is  equivalent  to  c  > 


or  t  a  +  y  . 


QED 


The  following  theorem  gives  a  useful  bound  on  ^(p)  for  small 
or  for  small  a. 


Theorem  2:  For  p  >  \  we  have  t^Cp)  >  (a+y)  v(p)  where 


>Kp)  =  2 


l°g(l^p)  ~  l°g(p) 


_1 _ 1 

1-P  ~  P 


) 


We  have  V(p)  <  1  for  p  >  and  lim  T(p)  =  1. 

P"2 


A  small  table  is 


given  below. 


1 


p 

Mp) 

.55 

.9933 

.60 

.9731 

.65 

.9389 

.70 

.8897 

.75 

.8240 

.7394 

.85 

.6319 

.90 

,4944 

For  a  =  y*  P  =  .80,  the  tables  of  Solomon  (1960)  give  tg(p)  -79  and 
the  bound  is  .7394. 

Proof :  We  repeat  the  argument  of  the  earlier  theorem  (proving 
tg(p)  < a  +  y)  but  with  the  inequalities  reversed. 

0 

To  show  -r—  F  (t)  >0,  it  suffices  to  show 
9p  P 


(*) 


h(s) 


+  log 


h(s) 


h(l-s) 


+  log 


h(l-s)l 


for  0  <  s  <  where  c  =  .  Observe  that  the  function  h(s)  is 

increasing  and  h(s)  +  h(l-s)  =  1.  Thus  lemma  5  part  (b)  says  that 
(*)  will  be  true  if 


h(0) 


+  log 


fh(0)l 


-  h(l) 


+  log 


h(l) 


Using  h(0)  =  1-p,  h(l)  =  p  and  a  little  algebra  this  becomes 


log<T^  "  lo^b 


14 


The  following  lemma  will  be  needed  in  the  proof  of  Theorem  3. 

Lemma  6 : 


I  ['tx(l-x)  +  a(l-2x)]x  (1-x)  e  dx  =  0  . 

'0 


Proof : 


IX 


\3i  , ,  .  a-1,,  a-1 

x  (1-x)  ]  =  a(l-2x)x  (1-x; 


so  that  integration  by  parts  gives 


[x3(l-x)'l]['teXtX]dx  =  - 


0 


[a(l-2x)x°l  ■'"(l-x)a  1]e^tXdx 


0 


Moving  both  integrals  to  the  same  side  of  the  equation  and  then  combining 
them  into  one  integral  completes  the  proof.  QED 


Theorem  3:  For  t  <  a  and  -i-  <  p  <  1, 


TZ  F  (t)  >  0  . 
3p  p 


Proof:  We  have  earlier  used  the  fact  (Lemma  1)  that  t—  F  (t)  has 

-  dp  p 


the  same  sign  as 


sp(t)  = 


,  .  a-1..  *a-l  Xtx. 

(p-x)x  (1-x)  e  dx 


where  X 


2p-l 


.  .  Therefore  it  suffices  to  show  that  S  (t)  >  0 
p(l-p)  p 


for  t  <  a.  This  can  be  achieved  by  noting  that  S  (t)  has  the  same 

P 

sign  as  e  S^(t)  where  $  is  any  smooth  function  depending  on  p 

and  t.  Now  (0)  is  positive  and  if  e^S^Ct)  is  increasing  in  t 

for  t  <  a,  then  sp(t)  must  be  positive  for  t  <  a.  So  it  suffices 

to  show  that  for  some  smooth  function  $  we  have  —  (e^S  (t))  >  0 

3t  p 

for  t  <  a. 

For  any  function  <p  we  have 


£  <*%<'»  ■ 


(p-x)  ( Xx+k)xa  1  (l-x;01  ^e  tXdx 


3: 


with  k  =  .  By  the  preceding  lemma  we  may  replace  (p-x)(Xx+k) 

in  the  above  integral  by 


(*) 


(p-x)(‘x+k)  +  c[\tx(l-x)  +  a(l-2x) ] 


without  changing  the  value  of  the  integral.  Here  c  is  any  quantity  whose 

value  does  not  depend  on  x  but  may  depend  on  t  and  p.  By  appropriate 

choice  of  c  and  k,  the  expression  (*)  can  be  made  into  a  constant, 

that  is,  an  expression  not  involving  x.  Choosing  c  =  -1/t  eliminates 
2 

the  quadratic  (x  )  term.  Next  we  choose  k  =  (2a/t)  -  X(l-p)  to 
eliminate  the  linear  term.  With  these  choices  the  constant  term  becomes 

pk  +  ca  *  (2p-l)(^-l)  >  0  for  t  <  a.  Thus  (e^S  (t))  >  0  for 

t  dtp 

A 

t  <  a  where  $  is  determined  by  *  (2a/t)  -  X(l-p).  QED 


3. 


POSITIVE  LINEAR  COMBINATIONS  OF  TWO  WE I BULL  VARIABLES 


In  this  section  it  is  shown  that  a  positive  linear  combination 
(-*1^1+7.^.,)  of  two  independent  identically  distributed  Weibull  random 
variables  Y^  has  a  distribution  function  which  is  Schur-concave  in  8 
for  sufficiently  large  argument.  Without  loss  of  generality  we  set 
8^+8-  =  1  and  we  set  the  scale  factor  of  the  Weibull  distribution  to 
one.  Then  ^  +  '  2^2  —  *s  ^chur-concave  in  1  if 

1 

(a)  t  _  (1  +  ~  when  6  >  1 

4m  ^ 

cr 

1  id 

(b)  t  >  (1  +  j?)3  2  ^  when  0  <  6  £  1  . 

Again  the  Schur-concavity  will  follow  if  the  derivative  with  respect 

to  p  of  the  distribution  function  P[pY^ +  (l-pJY^  t]  is  negative  in 
an  appropriate  region.  Theorem  4  developes  conditions  under  which  the 
derivative  is  negative  and  Corollary  1  interprets  them  to  provide  the 
result  above. 

Notation:  Define 

Fp(t)  =  P[pY1  +  (l-p)Y2  <  t]  . 

Theorem  4:  Let  Y^,Y2  independent  Weibull  random  variables  with 
6—1  6  1 

density  given  by  8xp  exp(-x  ).  If  +  <  p  <  1  and 


t"  >  <1+  77*  . 


then  ■=—  F  (t)  <  0.  Here 
3p  p 


for  B  >  1 


pB  +  (l-p)B  for  B  <  1 


Corollary  1:  Under  the  assumptions  of  Theorem  A,  we  have 


—  F  (t)  <  0  if 

3p  p 


(a)  t  _>  (1  +  -^-)B  when  B  _>  1 


1  IzA 

(b)  t  >_  (1  +  -~'>6  2  S  when  0  <  B  1  1  . 


Proof :  Part  (a)  follows  immediately  from  the  theorem.  It  is 

8  8 

necessary  to  find  an  upper  bound  for  p  +  (1-p)  independent  of  p 
for  part  (b).  Because  for  -^  <  p  <  1, 


{p6  +  (l-p)6} 


is  a  decreasing  function  of  p,  the  maximum  occurs  for  p  =  y,  i.e. 


{(■|)e  +  (^)6)  -  (y)B-1  -  21  6  . 

1  1^6 

Thus  tS  >_  (1  +  i)21-6,  i.e.  t  >  (1  +  -^-)B  2  6  .  QED 


The  following  lemma  is  used  in  the  proof  of  Theorem  4. 


Lemma  7 :  For  independent  Weibull  variables  with  density  given 

in  Theorem  4,  we  have 


exp 

-t  64>  (s )  1 

~  F  (t)  =  B2t26 

op  p  j 

~  .  8-1,,  v B" 1 
(l-2s)s  (1-s) 

0 

•  h(s) 

,  ,  ,2B+1 
h(s)  J 

ds 


C  g 

where  c (s)  =  s  +  (1-s)  and  h(s)  =  ps  +  (l-p)(l-s). 


Proof  of  Lemma:  Change  from  Cartesian  coordinates  to  r,s,  defined 

x 

bv  r  =  x+v,  s  =  — - — .  Then 
x+y 


fl  ft/h(s) 

F  (t)  =  !  [rf (rs)f (r(l-s)) ]drds 

p  J0  J0 


and 


4  F  (t)  =  f1 


9p  p 


u-2s)  -£-j  if(^) 


h(s) 


where  h(s)  =  ps  +  (l-p)(l-s)  and  f  is  the  Weibull  density.  QED 


Proof  of  Theorem  4:  Use  Lemma  7  and  let  g  denote  the  expression 
in  brackets  inside  the  integral  for  -r—  F  (t).  The  function 

dp  p 

o  i  8—1  1 

(l-2s)s  (1-s)  is  odd  about  the  value  s  =  y  .  Therefore,  to  prove 

d  1 

—  F  (t)  <  0,  it  suffices  to  show  that  g(s)  <  g(l-s)  for  0  <  s  <  -z  . 

dp  p  ^ 

Taking  logs  this  becomes 


X  |t-* 


8  8 

-  *-~Sl  +  (28+1)  log  h(s)  >  -  *  (1-~S^  +  (28  +  1)  log  h(l-s) 
h(s;  h(l-s)B 


or  equivalently 


c4> (s)  ,  ,  ,h(s)B  „  cq> (1-s)  rh(l-s)B 

h<s>*  108  ^  +  ^ 


8t 


6 


where  c  =  — .  To  obtain  this  last  expression  we  have  used  the 
fact  that  <f  (s)  =  4> ( 1— s ) .  By  Lemma  5(a)  (concerning  the  function 
+  log  x)  we  need  only  show 


h(s)g  h(l-s)6 
c<p(s)  cf(l-s)  * 


h(s)B  h(_l-s;B  < 
(11)  C?(s)  +  co(l-s)  -  2 


for  0  <  s  <  .  Note  that  (i)  is  true  because  h  is  strictly  increasing 

when  p  >  ^  and  (ii)  will  be  true  if 


Maximum 

i  '■ 

°ls±7 


[h(s)B  +  h(l-s)Bl  < 


s B  +  (l-s)S  1 


2c  . 


In  the  lemma  which  follows,  this  maximum  is  shown  to  equal  1  for  8  1 

6 


8  B  0 1 

and  to  equal  p  +  (1-p)  for  8  7  1.  Now  use  c  = 


28+1 


is  complete. 


and  the  proof 
QED 


Lemma  8:  Define 


'i'(s) 


h(s)B  +  h(l-s)^ 

8  B 

s  +  (1-s) 


for  6  >  1 


Maxi nun  '•‘(s)  = 
«  .  .  1 


1«  + 


[p  +  ( 1— p)  for  f?  <_  1  . 


Proof :  Remember  that  h(s)  -  ps  +  (l-p)(l-s)  and  —  <  p  <  1. 

Consider  first  E  '  1.  A  stationary  point  must  satisfy 

3 

log  v(s)  =  0  which  is  equivalent  to 


S-l  «_1  ft_1  ft_1 

£(2p-l) (h(s)"  -  h(l-s)  )  B(se  1  -  (l-sT  ‘) 


h(s)'  +  h(l-s)’ 


ft  ft 

sP  +  (1-sT 


I  Q  Q  O 

But  for  0  •'  s  <  4  it  is  easy  to  verify  that  h(s)  +  h(l-s)P>s  +  (1-s) 

and  h(s)3  1  -  hd-s)"  ^  <  sB  ^  -  (l-s)3  ^ .  Since  0  <  2p-j  <  1,  the 
left  hand  side  of  (*)  is  strictly  less  than  the  right  hand  side.  Thus 
has  no  stationary  point  in  the  open  interval  (0,4)  and  must 
achieve  its  maximum  at  0  or  . 


y(0)  =  p3  +  (l-p)S  >  1  =  if (4) 


Now  consider  °  >  1. 


~~  [h(s)B  +  h(l-s)S]  =  b(1"2s) [h (1-s) 3-1  -  h(s)6'1]  >  0 

dp  - 


13  1 

for  0  <  s  <  -z  .  Therefore  —  ui  (s)  >  0  for  0  <  s  <  -=■  .  The  subscript 

—  —  2  3p  p  —  —  —  2 

on  \p  indicates  the  value  of  p.  Since  ^(s)  increases  with  p  we 
have  if  (s)  <  1,(s)  =  1  =  41  (4)  as  desired.  QED 


4. 


POSITIVE  LINEAR  COMBINATION  OF  N  GAMMA  VARIABLES 


We  show  that  the  distribution  function  of  a  positive  linear  com¬ 
bination  of  n  independent  gamma  random  variables  is  Schur-concave  in 
the  coefficients  when  the  argument  of  the  distribution  function  is 
sufficiently  large.  This  implies  that  the  right  tail  probabilities 
are  Schur-convex.  The  "sufficiently  large"  bound  for  the  argument 
appears  to  be  too  large  when  compared  with  the  precise  results  of 
Section  2  for  n  =  2  and  can  probably  be  improved.  Corollary  2  shows 

— X 

that  for  independent  Y^  each  with  density  x  e  / F (ot )  and  positive 

n 

6.,  i  =  l,...,n,  P(  i_  6  .Y^  <_  t)  is  Schur-concave  in  0  for 
1  i=l  1 

n 

t  >_  (na  +  1)  1  9..  For  n  =  2  and  0..  =  p  and  0„  =  1-p  where 

i=l  1  1  1 

0  <  p  <  1  this  lower  bound  for  t  is  (2a+l)  which  is  too  large  by 
a  factor  of  2  according  to  the  results  of  Section  2. 

Theorem  5:  Let  ’  *  *^n  S3™03  random  variables 

with  density  xa  ^e  X/F(a)  and  9  =  (6  ,9  ...,9  )  where  9.  >  0 

l  L  n  i 

for  all  i.  Define 


n 

pf(e)  =  p[  l  9.Y  <  t]  . 
i=l  1  1 

As  a  function  of  9,  P ^  ( 0 )  is  symmetric  and  convex  (and  therefore  Schur 
convex)  inside  the  region  { 0 :  t  (na+l)Min  9^};  Pt(e)  is  symmetric  and 
concave  (and  therefore  Schur-concave)  inside  the  region  {9 :  t  >  (n<*+l)Max  9 


In  what  follows,  t  is  held  fixed  at  some  arbitrary  positive 
value.  The  basic  tool  in  the  proof  of  Theorem  5  is  the  following 
result  due  essentially  to  Marhshall  and  Proschan  (1965).  (See  Chapter 
11,  page  288  of  Marshall  and  Olkin  (1979).) 


Proposition:  If  X^,X0,...,X  are  exchangeable  and  g  is  a 

n 

continuous  convex  function,  then  ^  (a^  ,a^ , . .  .  .a^)  =  Eg(  )_  a^X^ )  is 

i=l 

symmetric  and  convex.  It  is  also  true  upon  replacing  convex  by  concave. 

n 

Proof  of  Theorem  5:  For  all  i  define  S.  =  Y./(  r  Y  ).  The 
-  r  i  kI1  k 

n 

vector  (S. , S.  , . . .  ,S  )  is  exchangeable  and  independent  of  ^  Y.  . 

i  1  n  .  “ ,  l 

1=1 

These  are  standard  properties. 

n 

Define  F(u)  =  P[  T  Y.  <  h] 
i=l  1 

•u  . 

i  na-1  -x,  . 

■  I  x  e  dx/F(na). 

•  0 


Then 


p[  l  O.Y.  1  t]  =  P[ (  )  0  ,S.)(  V  Y.)  <  t 
i=l  1  i=l  1  1  i=l  1 


EPfde.spdY.)  <  tlle.s.] 


=  EFCt/Ve.s.)  =  Eg^e.sp 


where  g(u)  =  F(t/u).  A  short  calculation  shows  that 


na  -t/u 

g"(u)  =  [  (na+1)  -  i]  ± —  . 

u  „  .  .  na+2 
(a)u 


Thus  g  is  convex  for  u  ■>  — and  concave  for  0  <  u  <  — — r  .  For 

—  na+1  —  na+1 


in  the  region  {?:t  (na+l)Min  0 . } , 


p.S.  >  -TT  • 
i  i  —  na+1 


p 

Thus  Eg(_0  S  )  is  symmetric  and  convex  by  the  proposition.  Similarly 


for  in  \-':t  (na+l)Max 


[  e.s.  < 
L  11  — 


and  therefore  EgC^^S^)  is  symmetric  and  concave.  QED 
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